A 2D Quantum Walk Simulation of Two-Particle Dynamics 
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Multi-dimensional quantum walks can exhibit highly non-trivial topological structure, providing 
a powerful tool for simulating quantum information and transport systems. We present a flexible 
implementation of a 2D optical quantum walk on a lattice, demonstrating a scalable quantum 
walk on a non-trivial graph structure. We realized a coherent quantum walk over 12 steps and 
169 positions using an optical fiber network. With our broad spectrum of quantum coins we were 
able to simulate the creation of entanglement in bipartite systems with conditioned interactions. 
Introducing dynamic control allowed for the investigation of effects such as strong non-linearities 
or two-particle scattering. Our results illustrate the potential of quantum walks as a route for 
simulating and understanding complex quantum systems. 



Quantum simulation constitutes a paradigm for de- 
veloping our understanding of quantum mechanical sys- 
tems. A current challenge is to find schemes, that can be 
readily implemented in the laboratory to provide insights 
into complex quantum phenomena. Quantum walks [I]- 
[3] serve as an ideal test-bed for studying the dynamics of 
such systems. Examples include understanding the role 
of entanglement and interactions between quantum parti- 
cles, the occurrence of localization effects [4], topological 
phases [5 , energy transport in photosynthesis [6] [7] , and 
the mimicking of the formation of molecule states [8]. 
While theoretical investigations already take advantage 
of complex graph structures in higher dimensions, exper- 
imental implementations are still limited by the required 
physical resources. 

All demonstrated quantum walks have so far been re- 
stricted to evolution in one dimension. They have been 
realized in a variety of architectures, including photonic 
[9HT2] and atomic [T3HT5] systems. Achieving increased 
dimensionality in a quantum walk [16 is of practical in- 
terest as many physical phenomena cannot be simulated 
with a single walker in a one-dimensional quantum walk, 
such as multi-particle entanglement and non-linear in- 
teractions. Furthermore, in quantum computation based 
on quantum walks [T7J [18] search algorithms exhibit a 
speed-up only in higher dimensional graphs [T9H22] . The 
first optical approaches to increasing the complexity of 
a linear quantum walk [23-25] showed that the dimen- 
sionality of the system is effectively expanded by using 
two walkers, keeping the graph one-dimensional. While 
adding additional walkers to the system is promising, in- 
troducing conditioned interactions and in particular con- 
trolled non-linear interactions at the single photon level 
is technologically very challenging. Interactions between 
walkers typically result in the appearance of entangle- 



ment, and have been shown to improve certain applica- 
tions, such as the graph isomorphism problem [26]. In 
the absence of such interactions, the two walkers remain 
effectively independent, which severely limits observable 
quantum features. 

We present a highly scalable implementation of an op- 
tical quantum walk on two spatial dimensions for quan- 
tum simulation, using frugal physical resources. One ma- 
jor advance of a two-dimensional system is the possibility 
to simulate a discrete evolution of two-particles including 
controlled interactions. In particular, one walker, in our 
case a coherent light pulse, on a 2D lattice is topologically 
equivalent to two-walkers acting on a one-dimensional 
graph. Thus, despite using an entirely classical light 
source, our experiment is able to demonstrate several 
archetypal two-particle quantum features. For our sim- 
ulations we exploit the similarity between coherent pro- 
cesses in quantum mechanics and classical optics [27] [28] , 
as it was used for example to demonstrate Grover's quan- 
tum search algorithm [29] . 

A quantum walk consists of a walker, such as a photon 
or an atom, which coherently propagates between dis- 
crete vertices on a graph. A walker is defined as a bipar- 
tite system consisting of position (x) and a quantum coin 
(c). The position value indicates at which vertex in the 
graph the walker resides, while the coin is an ancillary 
quantum state determining the direction of the walker 
at the next step. In a two-dimensional quantum walk 
the basis states of a walker are of the form \x u x 2 , c 1; c 2 ) 
describing its position x lj2 in spatial dimension 1 and 2 
and the corresponding two-sided coin parameters with 
c 1)2 = ±1. The evolution takes place in discrete steps, 
each of which has two stages, defined by coin (C) and 
step (S) operators. The coin operator coherently ma- 
nipulates the coin parameter, leaving the position un- 
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FIG. 1: (A) Experimental setup. Our photon source is a pulsed diode laser with pulse width 88ps, wavelength 
805nm and repetition rate 110kHz. The photons are initialized at position \x u x 2 ) = |0, 0) in horizontal polarization 
(corresponding to coin state |c l5 c 2 ) = | — 1,-1)). Once coupled into the setup through a low reflectivity beam 
splitter (BS, reflectivity 3%), their polarization state is manipulated with an EOM and a half- wave plate (HWP). 
The photonic wave packets are split by a polarizing beam splitter (PBS) and routed through single-mode fibres 
(SMF) of length 135m or 145m, implementing a temporal step in the x 2 direction. Additional HWPs and a second 
PBS perform a step in the x 1 direction based on the same principle. The split wave packet after the first step with 
equal splitting is indicated in the picture. At each step the photons have a probability of 12% (4%) in loops x 1 — 1 
(x 1 + 1) of being coupled out to a polarization and hence coin state resolving detection of the arrival time via four 
avalanche photodiodes (APDs). Including losses and detection efficiency, the probability of a photon continuing the 
walk after one step is 52% (12%) without (with) the EOM. (B) Projection of the spatial lattice onto a 
one-dimensional temporally encoded pulse chain for step one and two. Each step consists of a shift in both x 1 
direction, corresponding to a time difference of Ar x = 3.11ns, and x 2 direction with Ar 2 = 46.42ns. 



changed, whereas the step operator updates the position 
according to the new coin value. Explicitly, with a so- 
called Hadamard coin C H = H 2: a single step in the 
evolution is defined by the operators, 

H z \x^±l) (\x i: l)±\x i: -l))/V2 : Vi = l,2 
S\x 1: x 2 , c 1: c 2 ) y \x ± + c l5 x 2 + c 2 , c l5 c 2 ). (1) 

The evolution of the system proceeds by repeatedly ap- 
plying coin and step operators on the initial state 1^), 
resulting in |^) = (SC) n \ilj in ) after n steps. The step 
operator S hereby translates superpositions and entan- 
glement between the coin parameters directly to the spa- 
tial domain, imprinting signatures of quantum effects in 
the final probability distribution. 



We performed 2D quantum walks with photons ob- 
tained from attenuated laser pulses. The two internal 
coin states are represented by two polarization modes 
(horizontal and vertical) in two different spatial modes 
[App. 1], similar to the proposal in [3D]. Incident pho- 
tons follow, depending on their polarization, four differ- 
ent paths in a fiber network (Fig. 1A). The four paths 
correspond to the four different directions a walker can 
take in one step on a 2D lattice. Different path lengths 
in the circuit generate a temporally encoded state, where 
different position states are represented by discrete time- 
bins (Fig. IB). Each round trip in the setup implements 
a single step operation while the quantum coin opera- 
tion is performed with linear optical elements (half-wave 
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plates, HWP) [App. 1]. In order to adjust the coin oper- 
ator independently at each position we employed a fast- 
switching electro-optic modulator (EOM). A measure- 
ment with time-resolving single-photon counting mod- 
ules allowed for the reconstruction of the output photo- 
statistics [App. 2]. 

We have implemented two different kinds of quantum 
coins in our 2D quantum walks. First we investigated 
quantum walks driven only by separable coin operations, 
C = C x C 2 . Here the separability can directly be ob- 
served in the spatial spread over the lattice, when ini- 
tializing the walker in a separable state. As an example 
we measured a Hadamard walk with photons initially 
localized at position \x 1: x 2 ) = |0,0). The probability 
distribution showing at which position the photons were 
detected after ten steps (Fig. 2A±B) can be factorized 
into two independent distributions of one-dimensional 
quantum walks [IB] , stating no conceptual advantage of 
a 2D-quantum walk. However, two-dimensional quan- 
tum walks allow for much greater complexity using con- 
trolled operations. These operations condition the trans- 
formation of one coin state on the actual state of the 
other. Due to the induced quantum correlations one ob- 
tains a non-trivial evolution resulting in an inseparable 
final state. The probability distribution for a Hadamard 
walk with an additional controlling operation can be 
seen in Fig. 2C±D. We compare the ideal theoretical 
distribution with the measured photo-statistics via the 
similarity, S = ( V P ^( x n x 2 )P exp (x 1 , x 2 )) 2 , quantify- 

x ± ,x 2 

ing the equality of two classical probability distributions 
(5 = for a completely orthogonal distributions and 
S = 1 for identical distributions). For the Hadamard 
walk (Fig. 2A±B) we observe S = 0.957 ±0.003, and for 
the quantum walk with controlling gates (Fig. 2C+D) 
S = 0.903 ± 0.018 (after 10 steps, across 121 positions). 

Increasing the number of walkers in a quantum walk ef- 
fectively increases its dimensionality [23] . Specifically, for 
a given ID quantum walk with N positions and two walk- 
ers, there exists an isomorphic square lattice walk of size 
N 2 with one walker. By this topological analogy, a mea- 
sured spatial distribution from a 2D lattice with positions 
(x u x 2 ) can be interpreted as a coincidence measurement 
for two walkers at positions x 1 and x 2 propagating on the 
same linear graph. Hereby each combined coin operation 
of both particles, including controlled operations, has an 
equivalent coin operation in a 2D quantum walk. This 
allows us to interpret the 2D walk in Fig. 2C±D as a 
quantum walk with controlled two-particle operations, a 
system typically creating two-particle entanglement. The 
inseparability of the final probability distribution is then 
a direct signature of the simulated entanglement. 

In Fig. 2E we show a lower bound for the simulated en- 
tanglement between the two particles during the stepwise 
evolution with four different coin operations. We quan- 
tified the simulated entanglement via the von Neumann 



entropy E, assuming pure final states after the quantum 
walk [App. 4]. For this calculation the relative phases 
between the positions and coins were reconstructed from 
the obtained interference patterns, while phases between 
the four coin states were chosen to minimize the entan- 
glement value. Without conditioned operations the two 
particles evolve independently {E = 0), whereas an evo- 
lution including controlled operations reveals a probabil- 
ity distribution characterized by bipartite entanglement. 
We found that the interactions presented in Fig. 2C+D 
exhibit an entropy of at least E = 2.63 ±0.01 after 12 
steps, which is 56% of the maximal entropy (given by 
a maximally entangled state). The non-zero entropies 
obtained in the higher steps of the separable Hadamard 
walk are attributed to the high sensitivity of the entropy 
measure to small errors in the distribution for E w 0. 

The investigated interactions can be interpreted as 
long-distance interactions with the interaction strength 
being independent of the spatial distance of the parti- 
cles. This is a unique effect and highly non-trivial to 
demonstrate in actual two-particle quantum systems. 

Contrary to the position independent interactions is 
the evolution of two-particle quantum walks with short- 
range interactions, that is interactions occurring only 
when both particles occupy the same position. These 
interactions can be interpreted as two-particle scattering 
or non-linear interactions. When utilizing a 2D quantum 
walk to simulate two walkers, all vertices on the diagonal 
of the 2D-lattice correspond to both walkers occupying 
the same position. Hence, we can introduce non-linear 
interactions by modifying the coin operator on the di- 
agonal positions while keeping all other positions unaf- 
fected. As an example of a two-particle quantum walk 
with non-linear interactions (Fig. 3), the coin operator 
on the diagonal is in the form C nl = (H 1 <g> H 2 )C Zl where 
C z is a controlled phase operation implemented by a 
fast switching EOM. The chosen operation simulates a 
quantum scenario of particular interest - the creation 
of bound molecule states, predicted as a consequence 
of two-particle scattering [8 j. Evidently, the quantum 
walk is to a large extent confined to the main diagonal 
(%2P(x,x) = 0.317 ±0.006 as opposed to the Hadamard 

X 

walk J2P(x,x) = 0.242 ±0.001), a signature of the pres- 
ence of a bound molecule state. In general, using a coin 
invariant under particle exchange, bosonic or fermionic 
behavior can be simulated, depending on whether the ini- 
tial states are chosen to be symmetric or anti-symmetric 
with respect to particle permutations. With our initial 
state being invariant under particle exchange we simu- 
lated an effective Bose-Hubbard type non-linearity for 
two bosons [31] , We have demonstrated an efficient im- 
plementation of a two-dimensional quantum walk and 
proved the experimental feasibility to simulate a diversity 
of interesting multi-particle quantum effects. Our exper- 
iment overcomes the technical challenges of two-particle 
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FIG. 2: Measured and simulated probability distribution P(x 1 ,x 2 ) (traced over the coin space) after ten steps of a 
2D quantum walk with initial state |0,0, —1, —1). Theoretical (A) and measured (B) probability distribution of a 
2D Hadamard walk using the operation C H (Eq. 1). As only separable coin operations were performed (inset), the 
distribution is separable, given by a product of two one-dimensional distributions (gray). Theoretical (C) and 
measured (D) probability distribution of a 2D walk with controlled-Not X and controlled-phase operation Z, 
resulting in an unfactorizable distribution. Here c 2 is only transformed by XZ\ ±1) — >• ±| =F 1) if c x = —1 . The 
results (B) and (D) are obtained by detecting over 7 x 10 3 events and calibrated by the detection efficiencies of all 
four coin basis states. (E) Dynamic evolution of the Von Neumann entropy E generated by quantum walks (B) and 
(D) and quantum walks using controlled Hadamard coin operations (inset). The experimental values (dots) and 
theoretical predictions (dashed lines) mark a lower boundary for simulated two-particle entanglement. Statistical 
errors are smaller than the dot size. 
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FIG. 3: (A) Circuit representation of coin operations simulating non-linear interactions via 2D quantum walk. Only 
when the two virtual particles meet (x 1 = x 2 ) a controlled operation is applied. Theoretical (B) and measured (C) 
coincidence distribution P(x 1 ,x 2 ) (traced over the coin space) after seven steps of a simulated two-particle quantum 
walk with initial state |0, 0, —1, —1). The high probability that both particles are at the same position (diagonal) is a 
striking signature of bound states. The measured distribution is reconstructed by detecting over 8 x 10 3 events and 
has a similarity of S = 0.957 ± 0.013. Adding the EOM to the setup for dynamical control limits the step number to 
n = 7 due to the higher losses per step. Small imperfections of the EOM are included in the theoretical plot. 



experiments, while exhibiting very high similarity and 
scalability. Combined with the flexibility in the choice 
of input state, controlling the coin at each position inde- 
pendently allows for simulations of a broad spectrum of 
dynamic quantum systems under different physical con- 
ditions. 

Our experimental architecture can be generalized to 
more than two dimensions, with the addition of ex- 
tra loops and orbital angular momentum modes as coin 
states [32]. This opens a largely unexplored field of re- 
search, facilitating quantum simulation applications with 
multiple walkers, including bosonic and fermionic behav- 
ior, and non-linear interactions. It may be possible to 
study the effects of higher dimensional localization, graph 
percolations or utilize the network topology in conjunc- 
tion with single- or two-photon states. Additionally, a 
foreseeable future application for our system is the im- 
plementation of a quantum search algorithm. We demon- 
strated that, with a physical resource overhead, a classi- 
cal experiment can simulate many genuine quantum fea- 
tures. While our experiment is important for simulation 
applications, it is equally interesting for understanding 
fundamental physics at the border between classical and 
quantum coherence theory. 
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Appendix 1: Quantum gates with optical elements 

We realized the four internal coin states | ± 1, ±1) with 
the linear polarization states, horizontal \H) and vertical 
|V), and two spatial modes \a) and similar to the 
four spatial modes proposed in (27). The spatial modes 
correspond to the two input ports of the first polarizing 
beam splitter (Fig. A.l)). We encoded the states by 



|iT,a)-H-l,-l); |y,a)-H-l,+l); 
|ff,6)-^| + l,+l) ; \V,b) + 



(Al) 



To implement our quantum operations in the four- 
dimensional Hilbert space of the quantum coin we de- 
composed the U (4) unitary coin operation into products 
of multiple f7 (2) operations (31). Each U(2) transforma- 
tion is implemented either by half- wave plates (HWP) or 
an electro-optic modulator (EOM). 

In the basis of the four coin states the transformations 
are given by 



HWPx 



( cos(2(9 1 ) 
sin^) 


V 



sin^) 0\ 

-cos(2(9 1 ) 

1 

1/ 



HWP o — 
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HWP 3 



/ 1 

10 

cos(2<9 2 ) sin(2(9 2 ) 

\00 sin(20 2 ) -cos(2(9 2 ) / 

/cos(2<9 3 ) sin(2<9 3 ) \ 

10 

1 

\ sin(2(9 3 ) -cos(2f9 3 ) / 

/ 1 \ 

-cos(2<9 4 ) sin(2(9 4 ) 

sin(2<9 4 ) cos(2<9 4 ) 

\ 1 / 

/ e .*(*i,* 2 ) \ 

10 

10 

\ 1 / 



(A2) 



with Q i being the angle of HWP^, i = {1,..,4}, relative 
to its optical axis and <fr(x u x 2 ) a tunable phase. 

Depending on the position of the four HWPs in the 
setup (Fig. A.l A) they operate on different coin state 
pairs, due to the spatial switch via polarizing beam split- 
ters (PBS). The Hadamard coin (Fig. 2B+E) was ob- 
tained for the configurations { = 7r/8,Vz, while the coin 
used in Fig. 2D+E was given by 1 = — 7r/8, 2 _ 4 = ir/8. 
To implement controlled-Hadamard gates either HWP 2 3 
or HWP 2 4 were aligned to their optical axis (6 = 0), 
while the remaining plates were set to = tt/S (Fig. 2E). 



The transformation of the EOM with <j)(x 1 ,x 2 ) = tt cor- 
responds to a controlled-Z operation. 

Additional static phase factors changing the relative 
phase between the four coin states can occur during the 
propagation through the setup. However, these phases do 
not influence the final propability distribution due to the 
property of the coin operators that can be implemented 
with the used optical elements. Given precise phase con- 
trol, two additional HWPs and a PBS would allow the 
implementation of arbitrary U (4) coin operators. 



Appendix 2: Quantum walk implementation via 
time- multiplexing 

Our experiment simulates a 2D quantum walk on a 
regular square lattice, which means that a walker can 
move in four possible directions from a given site. The 
direction of the movement is determined by the current 
coin state of the walker. To implement the quantum 
walk in a 2D topology we use the time-multiplexing tech- 
nique (9). This method maps each individual position of 
the 2D graph on the one-dimensional time line. In con- 
trast to determining the direction of the following step 
in space, the coin state defines a fixed time delay in the 
time-multiplexed system. 

Fig. A.l shows the first step of a quantum walk in our 
system and the corresponding mapping of the temporal 
spread on the 2D lattice. The propagation of the initial 
pulse representing the walker through the setup is split 
into five parts: 

(1) We initialize the photonic input state in the coin 
state \H,a) with the starting time corresponding 
to position \x x ,x 2 ) = |0,0). 

(2) After launching the pulse into the setup we per- 
form the coin operation C 2 = C HW p 2 C'hwp 1 C7eom, 
as defined in Eq. A2. 

(3) Split by a polarizing beam splitter (PBS) the pulse 
travels through two fibers of different lengths. The 
resulting time separation Ar 2 can be interpreted as 
a step in vertical (x 2 ) direction. We hereby define 
the transformation x 2 — >> x 2 — 1 for pulses passing 
the longer fiber, while the pulses in the shorter arm 
acquire the transformation x 2 — » x 2 + 1. Additional 
retardation plates hereby compensate for unwanted 
polarization effects in the fibers. 

(4) Subsequently HWP 3 and HWP 4 implement the 
coin operation C 1 = C HWP3 C HWP4 (Eq. A2), deter- 
mining the coin state for the horizontal (x ± ) direc- 
tion. 

(5) The step operator in the horizontal (x x ) direction is 
performed with a second PBS and a temporal shift 
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FIG. A.l: (A) Schematic setup including five time stages during the first step of a Hadamard walk. The laser pulse 
is initialized at (1) in horizontal polarization (arrow). After the coin operation C 2 (2) the vertical step (3) is 
implemented via two single-mode fibers (SMF) of different lengths. Subsequently C 1 (4) performs the coin operation 
for the horizontal direction, followed by the step operation implemented in free space (5). We compensate for 
different losses in arms a and b via different splitting ratios of the according BS. The corresponding spatial spread 
on a 2D lattice is shown for part (1), (3) and (5). (B) Mapping of the timing information of detected events onto 
spatial coordinates. The minimal time for a round trip T min is 676ns. For more details see text. 
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At 15 obtained by traveling in two free space paths 
of different lengths. 

As a result one step in our quantum walk setup combines 
both a shift in vertical and horizontal direction, restrict- 
ing the translation on a 2D grid to the diagonal neighbors 
(Fig. IB). Starting in the origin, this leads to a spread 
over maximally (n + l) 2 positions after n steps. 

After each step a time resolving measurement gives 
information about the location of the photon. The de- 
tection times for the first and second step of a Hadamard 
walk are shown in Fig. A. IB. At the first detection the 
photons can arrive at four different times, which is either 
the minimal time T min (corresponding to position state 
|1,1», T min + At, stating position T min + Ar 2 

(| - 1, 1)) or T min + Ar ± + Ar 2 (| - 1, -1)). At the second 
step on the other hand, multiple pulses can arrive simul- 
taneously at the detectors. Independent of the coin state 
and, hence, spatial and polarization modes of the pho- 
tons, all coinciding detections correspond to the same po- 
sition state \x 1 ,x 2 ). In the following step wave packets in 
the same time bin and spatial mode interfere at HWP 12 
(C 2 ), while wave packets in different spatial modes can 
interfere at HWP 3 , 4 (Ci). 

Appendix 3: Limits and imperfections 

During the time evolution, the area of the grid cov- 
ered by the quantum walker grows quadratically with 
the number of steps. The use of the time- multiplexing 
technique guarantees that the number of elements stays 
constant independent of the size of the simulated grid. In 
the experimental implementation it is only the lengths of 
the optical paths that needs to be adjusted to the maxi- 
mum number of steps that are to be realized. In addition, 
the performance of our time-multiplexed setup is limited 
only by imperfections of the optical components result- 
ing in errors, decoherence and losses. In the following we 
want to discuss each point individually. 

A sharp limitation for the maximal step number is 
given by the design of the experiment. If the minimal 
time for one round trip T min is shorter than the temporal 
expansion of all positions in a single step, temporal over- 
laps between different steps can occur. We choose our 
experimental parameters T min and Ar 1>2 without EOM 
such that At 2 > 13 • Ar x and T min > 13 • Ar 2 to pre- 
vent temporal overlaps of different positions for the first 
12 steps. An occurring additional delay induced by the 
modulator changed the conditions without inducing un- 
wanted temporal overlaps. However, by simply changing 
the fiber lengths or the path differences the step number 
can easily be increased. 

The most significant source of systematic errors in the 
setup is the EOM. Due to the architecture of the mod- 
ulator the applied phase is not only affecting the hori- 



zontal polarization component, as shown in Eq. A2, but 
also the vertical component with a factor 1/3.5. This 
decreases the achievable similarities for quantum walks 
where controlled- Z coin operations are used. Addition- 
ally, the wave front of pulses passing the modulator are 
distorted differently for both polarizations, which influ- 
ences the occurring interferences. Both effects are in- 
cluded in the theory presented in Fig. 3B. A replacement 
of the EOM with an optimized modulator would improve 
the achievable similarity at higher number of steps. 

Decoherence in the time-multiplexed setup can occur 
if mechanical vibrations of the optical elements influence 
the interference properties. Typically mirrors vibrate 
with a frequency below 500Hz, corresponding to a time 
scale of 2ms. The duration of twelve steps in the current 
setup is less than 10/i«s, a factor of 200 faster compared to 
mirror vibrations. This suggests that decoherence effects 
will not influence the time-multiplexed quantum walk up 
to at least 100 steps. 

At the present stage the main factor limiting the scal- 
ability is given by the losses per step. These are induced 
by the probabilistic detection method and losses at opti- 
cal elements. To counter the effect of losses one can either 
start with an increased intensity or use optical amplifiers, 
as shown in (32). While the first approach requires an ac- 
tive protection of the single-photon detectors, the second 
prohibits the use of the experiment with single-photon 
sources. A third method to reduce the losses is a change 
from probabilistic to a deterministic coupling mechanism 
with additional polarization modulators. This technique 
combined with a change-over to a low-loss wavelength 
regime (1550nm), makes the setup interesting for single 
photon input states. Using one of the described methods 
to circumvent the losses can increase the number of steps 
significantly. 

Appendix 4: Entanglement 

To quantify the two-particle entanglement simu- 
lated in the system we assumed that the quan- 
tum walk evolution results in a pure state \i/j n ) = 
S o, x1jX2jC1jC2 \x u x 2 ,c u c 2 ) after n steps, with the 

x\ ,x 2 ,c 1 ,c 2 

complex parameters a X1)X2iC1]C2 G C. The assumption is 
based on the fact that the system does not show signs of 
decoherence for any of the coin operators, confirmed by 
the high values of the measured similarities. 

The von Neumann entropy E, which quantifies the en- 
tanglement (33) is given by 

^(Pi) = -^A i log 2 A i , (A3) 

i 

with the eigenvalues of the reduced density matrix 
= Tr 2 (p) = Tr 2 (|'0 n )('0 n |), given by the trace over one 
subsystem. For details see (23). 
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By individually measuring the arrival probability in 
each coin and position state we obtain the absolute 
squared of the parameters \a xliX2tCliC2 1 2 , hence no direct 
extraction of the phase information is possible. However, 
we can obtain information about the relative phases be- 
tween position and coin states from the final interference 
pattern. As a result we can reconstruct the phase dis- 
tribution with the help of the theoretical model up to 
three undetermined relative phases between the four dif- 
ferent coin states and a global phase factor. By choosing 
the phases inducing the minimal entropy, we are able to 
give a lower bound for the simulated entanglement in the 
experiment. 
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